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Abstract. We explore the role of mass loss and convective core overshoot in the 
evolution of Classical Cepheids. Stellar evolution models are computed with a recipe 
for pulsation-driven mass loss and it is found that mass loss alone is unable to account 
for the long-standing Cepheid mass discrepancy. However, the combination of mass 
loss and moderate convective core overshooting does provide a solution, bringing the 
amount of convective core overshooting in Cepheids closer to that found for other stars. 



1. Introduction 

Classical Cepheids are powerful laboratories for understanding stellar astrophysics, yet 
an important unanswered question is what are the masses of Cepheids. Cepheid masses 
are determined using multiple different methods, including stellar evolution models 
matching measured effective temperatures and luminosities, stellar pulsation models 
matching pulsation periods and amplitudes, and from observations of binaries where 
one component is a Cepheid. However, these three methods do not agree such that 
masses based on stellar pulsation models tend to be smaller than masses based on stellar 
evolution m odels. Furthermor e, dynamic masses are consistent with stellar pulsation 
masses (e.g. lEvansetail 2008). 



This mass diff erence, called the Cepheid Mass Discrepancy, is a long-standing 
problem (ICoxlll980 ) that is important not just for understanding Cepheids themselves 
but how stars evolve in general. Because dynamic masses agree with the pulsation 
masses but not the stellar evolution masses then there must be physics m issing from the 



stella r evolution m odels leading to Cepheid masses being overestimated (Kell er & Wood 



20061 ; iKellerl 12008). By finding the underlying source of the mass discrepancy, we can 



constrain the evolution of massive stars in general, leading to a better understanding of 
their later phases, such as asymptotic giant branch stars, white dwarfs and supernovae. 
Historically, the Cepheid Mass Discrepancy was foun d to be abo ut 40% in gen- 



eral an d about a factor of two for beat Cepheids (ICoxlll980h . However. [Moskali k et al. 



1992b reduced the mass discrepancy substantially when the OPAL opacities were in- 



cluded in stellar evolution models. This result was a tremendous step forward in re- 
solving the mass discrepancy. However, ICaputo et all (I2005h and IKellerl (l2008h found 



that the mass discrepanc y for Galactic Cepheid s is still present, with typical values of 
10 - 20%. Furthermore. Keller & Woodl (l2006h modeled stellar pulsation of Galactic, 
and Large and Small Magellenic Cloud Cepheids and found that the mass discrepancy 
is a function of metallicity, increasing with decreasing metallicity. The purpose of this 
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work is to explore some of the physics used in stellar evolution calculations and test the 
impact of different input physics on the Cepheid mass discrepancy. 



2. Possible Causes of the Mass Discrepancy 

There have be en a number of possible solutions suggested to resolve the mass dis- 
crepancy (see lBono et al.ll200q. for more details). The four most probable causes are: 
miss ing opacity, rotatio nal mixing, convective core overshooting (CCO), and mass loss. 
Bono et all (2006) argued that for opacity changes to account for the discrepancy, 



they would need to differ by a factor of two. The difference between the Opacity Project 



(Badn ell et all 2005) and previous opacities is at most 10%, suggesting that an opacity 
revision of a factor of two appears unlikely. Therefore new opacities would only have 
a small impact on the mass discrepancy, and thus missing opacity is not a plausible 
solution. 

Rotational mixing is a second possible solution. A star that rotates rapidly dur- 
ing its main sequence evolution will mix hydrogen into the core. Thus its post main 
sequence helium core is more massive than if the same star evolved with negligible ro- 
tation. The more massive helium core would cause the star to be more luminous when 
the model crosses the Cepheid Instability Strip, and thus would predict a smaller mass 
for a Cepheid model with the same luminosity as a Cepheid model that does not include 
rotational mixing. Furthermore, rotational mixing will change the surface abundances 
of Cepheids as well as the structure of the post main sequence blue lo op evolution where 
stars undergo their second and third crossings of the Instability Strip (Maeder & Mevnet 
2001b . Rotational mixing is a possible solution to the mass discrepancy that should be 
explored fur ther, h owever, it is beyond the scope of this work. 

iKellerl d2008l and references therein) suggested convective core overshoot (CCO) 
as yet another possible solution. Convective core overshooting during the main-sequence 
evolution mixes hydrogen into the nuclear-burning region. Just like rotational mix- 
ing, this will lead to a more massive helium core and a more luminous Cepheid. The 
amount of CCO in stellar evolution models is often determined from the parameteriza- 
tion A c = a c Hp, where Hp is the pressure scale height, a c is a free para meter, and A r is 



the dis tanc e above the c onvective region that overshooting penetrates. Keller & Wood 



(120061) and IKellerl (2008) xmnd that to solve the mass discrepancy a c must equal 0.5 



1. However, this range of a c is larger than that found for other stars. 

Mass Loss occurs throughout the entire evolution of a star and hence suggests an- 
other possible solution to the mass discrepancy. However, mass loss during the main 
sequence and red giant stage of evolutio n of intermediate mass stars appears to negli- 



gibly contribute to the ma ss discrepancy (Lanz & Cat alal 19921 : IWillsoij|2000l) with the 
exception of the Be stars. iNeilson & Lestei (2008, 2009) developed a prescription for 



pulsation-enhanced mass loss during the Cepheid stage of evolution and found that the 
predicted Cepheid mass-loss rates are large enough to contribute to the mass discrep- 
ancy. However, it is uncertain if mass loss completely solves the mass discrepancy. 
Further evidence of Cepheid mass loss includes ob servations of asymmetry of Ha line 
profiles of Galactic Cepheids dNardetto et al.ll2008l) a s well as from mode lling the in- 
frared excess of Large Magellanic Cloud Cepheids dNeilson et al.ll2010l) . However, 



iKelleil (120081) noted that mass loss needs to be as efficient for low-mass Cepheids as for 
large-mass Cepheids, which would be seem unlikely. 
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3. Methodogy 

We test if mass loss during the Cepheid stage of stellar evolution is a plausible solution 
to the Cepheid mass discrepancy, either on its own or when coupled with various values 
for convective core overshoot, a c . We compute grids o f stellar evolution m odels using a 
one-dimensional hydrodynamic stellar evolution code (He ger et all 2000) for stars with 
initial mass M = 4,6, 8, and 9 M Q . These models are evolved with four different sets 
of input parameters: 

1. no pulsation-driven mass loss and no CCO, 



2. pulsation-driven mass loss, using the prescription from lNeilson & Lesterl (|2008) 
and no CCO, 

3. pulsation-driven mass loss and CCO, with a c = 0.1, and 

4. pulsation-driven m ass loss and CCO , with a c = 0.335, the value found by Brott 
et al. (in prep) (see VinketaL 2010, for details). 



We compute the contribution of mass loss to the mass discrepancy (AM/M) by subtract- 
ing the initial mass of the stellar model from the mass at the end of Cepheid evolution, 
and divide by the initial mass. We calculate the contribution of CCO towards the mass 
discrepancy to be that a c = 0.1 leads to AM/M = 2.5%, a c = 0.2 is AM/M = 5% and 
a c = 0.335 is AM/M = 8.375% because of the change of luminosity due to increasing 
the amount of CCO. The pulsation mass is computed from the luminosity, period, and 
other pulsation quantities; therefore we need to include the contribution of CCO in our 
computed mass difference. We adopt the mass of the stellar mode l at the end of Cepheid 
evolution as equivalent to the pulsation mass that was found by iKellerl (120081) and the 
initial stellar m ass with zero CCO i s equivalent to the evolution mass based on the def- 
inition used bv lCaputo et aD(l2005h andlKeller (2008). Therefore, this measurement is 



a prediction of the maximum contribution of the chosen input physics on the Cepheid 
mass discrepancy. Note that we assume pulsation-driven mass loss occurs only on the 
Cepheid Instability Strip. 



4. Results 

The stellar evolution tracks for each case and mass are shown in Figure Q] and the 
computed contribution to the mass discrepancy is shown in Tabled] For the a c = 0.335 
case, the 6, 8, and 9 M Q stars lose enough mass during th e Cepheid stag e of evolution 
to be consistent with the mass discrepancy determined by iKeilerl (120081) . Every stellar 
model undergoes blue loop evolution, where the width of the blue loop is determined 
by effective temperature range spanned, increases with initial mass implying that the 
stars undergo the Cepheid stage of evolution for each case. We do not model the radial 
pulsations but assume that the star is undergoing Cepheid evolution if the effective 
temperature and luminosity of the stellar model is inside the Cepheid Instability Strip. 



Once a stellar model enters the Instability Strip, we use the iNeilson & Lester 



(2008) prescription for mass loss where radial pulsation in the star generates shocks 
that enhance the wind. In this model, the mass-loss rate is determined by the stellar 
mass, luminosity, effective temperature, an d pulsation p eriod, where the period is com- 



puted from a Period-Mass-Radius relation (IGierenlll989r) . This relation has a predicted 
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Figure 1. Stellar evolutionary tracks for models with initial mass M = 
4,6,8,9 M Q and for the four scenarios discussed in the text. (Upper Left) no 
pulsation-driven mass loss nor CCO, (Upper Right) with pulsation-driven mass loss 
and zero CCO, (Lower Left) with pulsation-driven mass loss and CCO, a c =0.1, and 
(Lower Right) with pulsation-driven mass loss and CCO, a c = 0.335. The dashed 
lines are the borders of the Cepheid Instability Strip. 



period error of approximately 25% whi ch may lead to a mass-l oss rate uncertainty of 
approximately a few hundred percent. iNeilson & Lester! ( 20081) found that mass-loss 
rates of Galactic Cepheids ranges from 10 -11 - 1 0" 7 M^yr' 1 and th at mass-loss rates 
are potentially increasing with lower metallicity dNeilson & Lesterll2009l) . We apply 
this prescription to stellar evolution models and explore how mass loss affects blue 
loop evolution. 

Consider first the 8 M Q stellar evolution models that are shown in Figure |2] The 
tracks for models with and without Cepheid mass loss and zero CCO have blue loops 
spanning identical ranges of effective temperature and luminosity suggesting that mass 
loss alone has little affect on the evolution of the star. When CCO is considered then 
the location of the blue loop shifts to larger luminosities. For example, the blue loop of 
the stellar evolution model with a c = 0.1 has a luminosity that is about 0.1 dex larger 
than the models without CCO while the model with a c = 0.335 is about 0.25 dex more 
luminous than the models without CCO. 

However, it is not clear how mass loss changes with different amounts of CCO 
and thus how mass loss affects structure of the blue loop. For instance, the mass- 
loss rate depends on the ratio of luminosity and mass, L/M, which increases with a c . 
On the other hand, the predicted mass-loss rate increases with decreasing pulsation 
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Table 1 . The contribution of mass loss and CCO to the Cepheid Mass Discrepancy. 



Initial Mass (Af ) 


Case 1 


Case 2 


Case 3 


Case 4 


4 


0% 


6.25% 


6.00% 


8.72% 


6 


0% 


3.00% 


8.30% 


18.54% 


8 


0% 


1.75% 


4.50% 


16.24% 


9 


0% 


1.67% 


4.40% 


15.00% 



period which in turn increases with stellar radius. The radius increases with luminosity 
and thus with increasing a c . Thus a change of CCO leads to changes in fundamental 
parameters that work against each other in determining the mass-loss rate. In Figure |2j 
the only instance where mass loss appears to affect the evolution of the star is for the 
case where a c = 0.335. 

We explore this further in Figure [3) where we compare the predicted mass-loss 
rates for the 8 M Q stellar evolution models for three cases where pulsation-driven mass 
loss is assumed for the second and third crossings of the Instability Strip; the first 
crossing is too short to contribute a significant mass change in the 8 models. During 
the second crossing, the timescale is about twice as long for a c = 0.335 as opposed to 
a c = and 0.1, yet there are little difference in mass-loss rates. However, during the 
third crossing, the a c = 0.335 model has a crossing time that is almost an order-of- 
magnitude longer than the cases with less CCO and has an average mass-loss rate that 
is a factor of 2-3 larger than for the other cases. 

The predicted average mass-loss rate tends to increase with a c when the star is 
evolving along the Instability Strip. For increasing CCO, then the ratio of the luminosity 
and mass, L/M, increases for Cepheids. Hence the amount of radiation pressure acting 
on a wind is larger and plays a bigger role in drivin g the wind than the chan ge of 
pulsation period. This is consistent with the results of iNeilson & Lester! (|2009). The 
change in timescale that the 8 M Q stars spend on the Cepheid Instability Strip is also 
not directly due to the amount of CCO but on the mass-loss rate and its dependence 
on CCO. We hypothesize that the widt h of the blue loop decrease s if the mass-loss rate 
is large, similar to the result found by iBrunish & Willsonl d!987h . thus "trapping" the 
star in the Cepheid Instability Strip. During the second crossing, a Cepheid's effective 
temperature is increasing and its radius is decreasing with time. However, mass loss 
acts against the contraction of the stellar envelope, and "puffs up" the envelope. The 
more inflated envelope prevents the effective temperature from increasing at the same 
rate as the cases with smaller CCO and mass loss. 

We find that the 9 M Q models appear to behave similarly as the 8 M Q models. The 
width of the blue loop of the 9 M Q models is smaller for the a c = 0.335 case and is 
also trapped in the Instability Strip. The stellar evolution tracks appear unchanged for 
smaller values of a c . The 4 and 6 M Q stellar models behave differently. The models 
with a c = and zero pulsation-driven mass loss appear to undergo multiple blue loops 
and when pulsation mass loss is included the multiple loops disappear and the one blue 
loop has a smaller width for both models. The 4 and 6 M Q blue loops are more sensitive 
to mass loss relative to the 8 and 9 M Q . 

The 4 and 6 M© model blue loops are also sensitive to CCO. The blue loops for 
the a c = 0.1 models span a larger range of effective temperatures than the models with 
0. However, the a c = 0.335 models have a smaller blue loop, so much so that 
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Figure 2. Comparison of the 8 M Q stellar evolution tracks for the four scenarios 
explored in this work. The location and structure of the blue loops change due to the 
different input parameters. 



the 4 M Q model blue loop does not enter the Cepheid Instability Strip. This behavior 
is associated with the minimum mass at which stellar evolution models can form a 
blue loop. The minimum mass increases with increasing a c , such that at a c = 0.335 
the 4 M Q model is approximately the minimum mass. The 6 M Q model blue loop 
appears truncated at a c = 0.335 due to the combination of mass loss and convective 
core overshooting just like the 8 and 9 M Q models. 



5. Summary 



The results of this work suggest that the Cepheid ma ss discrepancy can be r esolved 
by the combination of pulsation-enhanced mass loss dNeilson & Lesterl l2008h on the 
Cepheid Instability Strip and convective core overshoot in the main sequence progen- 
itors. On its own, CCO is not a plausible solution because of m easurements of the 
amount of CCO in other stars. While this w ould require a c _ > 0.5 dKellerl I2008I) , from 
eclips ing binaries, the value of a c = 0.2-0. 4 dClausen et al. 2010I: Sandberg Lacy et all 
201fJ) . from ft Cephei sta rs a r = 0.28 + 0.1 (Lovekin & Goupill l2010h while from early 



B-type stars a c = 0.335 (jHu nteret al.ll2008l : IVink et al.ll2010h . Furthermore, large val- 
ues of a r may act to suppress the blue loop evolution in the most massive Cepheids. 
iKippenhahn & Weige rt ( 1990) noted that whether stars go through blue loop evolution 
depends on the mass and radius of the stellar core, the mass and change of hydrogen 
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Figure 3. The pulsation-driven mass-loss rates predicted for the 8 M models 
for the second and third crossings of the Cepheid Instability Strip and for the three 
scenarios that include pulsation-driven mass loss. The pulsation-driven mass-loss 
rate for the 8 M© during its first crossing is not shown because the change of mass is 
too small to contribute significantly. 



abundance of the hydrogen-burning shell. A more massive stellar helium core, such as 
that produced when CCO is included, acts to quench blue loop evolution. 

Mass loss, on its own, also is not a plausible solution because the Cepheid crossing 
timescales are too short, and/or the mass-loss rates are too small. However, the com- 
bination of CCO, with ct c = 0.2 - 0.4, and pulsation-driven mass loss together provide 
a plausible solution. This work is preliminary and we intend to explore a larger mass 
range up to M = 15 M Q as well as Small and Large Magellanic Clo ud metallicities 
to test the metallicity dependence of the mass discrepancy found by iKeller & Wood! 
(120061) . 
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